Dynamics of entanglement between two atomic samples with spontaneous scattering 
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We investigate the effects of spontaneous scattering on the evolution of entanglement of two atomic 
samples, probed by phase shift measurements on optical beams interacting with both samples. 
We develop a formalism of conditional quantum evolutions and present a wave function analysis 
implemented in numerical simulations of the state vector dynamics. This method allows to track 
the evolution of entanglement and to compare it with the predictions obtained when spontaneous 
scattering is neglected. We provide numerical evidence that the interferometric scheme to entangle 
atomic samples is only marginally affected by the presence of spontaneous scattering, and should 
thus be robust even in more realistic situations. 
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I. INTRODUCTION 

In recent years much interest has been devoted to the 
study of quantum entanglement^ one of the most pro- 
found consequences of quantum mechanics and by now 
thought of as a fundamental resource of Nature, of com- 
parable importance to energy, information, entropy, or 
any fundamental resource [l| . Entanglement plays a cru- 
cial role in many fundamental aspects of quantum me- 
chanics, such as the quantum theory of measurement, 
decoherence, and quantum nonlocality. Furthermore, en- 
tanglement is one of the key ingredients in quantum in- 
formation and quantum computation theory, and their 
experimental implementations 

Many theoretical and experimental efforts have been 
recently devoted to create and study entangled states of 
material particles by exploiting photon-atom interactions 
[^1^. The machinery of quantum non-demolition (QND) 
measurements provides an important tool toward this 
goal. In particular it is possible to probe nondestruc- 
tively the atomic state population of two stable states 
of an atomic sample by phase-shift measurements on a 
field of radiation non-resonantly coupled with one of the 
atomic state 0, IE IE S 13 • The consequent modification 
of the state vector or density matrix of the sample can 
be used to entangle pairs of samples by performing joint 
phase-shift measurements on light propagating through 
both samples that provide information about the total 
occupancies of various states |E]- Recently, entangled 
states of two macroscopic samples of atoms were realized 
by QND measurements of spin noise ,10i|. 

So far, most of the theoretical work has focused on 
the "static" entanglement properties, but very few re- 
sults have been obtained on the dynamical properties of 
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the entanglement of these systems, i.e., on how to real- 
ize the most efficient entangling dynamics by using the 
interaction processes allowed by the physical set-ups. 

In order to better understand the dynamics of entan- 
glement in coupled systems of matter and radiation, a 
quantitative wave function analysis has been recently 
introduced to determine the entanglement created 
by measurements of total population on separate atomic 
samples by means of optical phase shifts. The authors 
of Ref. ,11] considered a photon-atom interaction scheme 
in which dissipative effects, as spontaneous scattering, 
were completely neglected. However, a proper inclusion 
of spontaneous scattering would provide a first step to- 
ward an understanding of the decoherence effects affect- 
ing the system, and, in view of achieving a practical con- 
trol of the entangling protocol, it must be considered in 
order to have a more realistic physical description of the 
process. 

In this paper we reformulate the interaction model pre- 
sented in Ref. so as to include spontaneous scatter- 
ing in the wave function analysis, and provide a com- 
parison of the results obtained from the two schemes of 
photon-atom interaction, with and without spontaneous 
scattering. A byproduct of our study, that is not re- 
stricted to the specific problem faced in the present work, 
is the formulation of a theory of conditional quantum evo- 
lution for generic photon-atom scattering in systems with 
many atoms. 

The paper is organized as follows. In Sec^lwe present 
the interferometric set-up and the detection scheme used 
to measure the field phase shift and to detect the scat- 
tered photons. Moreover, we introduce the analytical 
model to determine how the wave function of the atomic 
samples is modified by the photo-detection. We will dis- 
tinguish between the case in which the photon is scat- 
tered in the same or in a different mode with respect 
to the initial one. In Sec IIIII we analytically show how 
the two atomic clouds get entangled due to the photo- 
detection and give the algorithm to implement the phys- 
ical model in a numerical simulation. The results of two 
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kinds of simulations are presented in Sec. IIVI In the 
first one, we simulate the consecutive measurement of two 
combinations of spin variables, showing how the entan- 
glement of the two atomic samples evolve as the photo- 
dctcction occurs. In the second one, we present simula- 
tions where the atomic samples are subject to continuous 
spin rotations during measurement, resulting in a differ- 
ent evolution of the entanglement. For both schemes of 
simulations we then compare the case in which the spon- 
taneous scattering is considered with the case in which it 
is neglected. Finally, in Sec. we draw our conclusions. 

II. THE INTERFEROMETRIC SCHEME AND 
THE WAVE FUNCTION UPDATING MODEL 

Let us consider two atomic samples placed in one of 
the arms of an intcrfcromctric set-up used to measure the 
phase shift of the entering electromagnetic field (Fig.^. 

A photon impinging the 50-50 beam splitter 1 is "sep- 
arated" in two components, denoted by F, which follows 
the upper free path, and /, which can interact with the 
two atomic samples in the lower arm. ap and a/ denote 
the annihilation operator modes of the two field compo- 
nents. 




FIG. 1: Atoms occupying the internal state \a) in the two 
samples interact with the light field which is incident from 
the left in the figure. The phase shift of the light field due 
to interaction with these atoms is registered by the different 
photo-currents in the two detectors. 

Each sample is composed of N atoms, whose level 
structure consists of two stable states, \a) and and one 
excited state, |c) (Fig.EJ. States \a) and |c) are coupled 
off-resonantly by the electromagnetic field whose annihi- 
lation operator is aj. Thus, a photon passing through the 
two atomic samples can be either absorbed and sponta- 
neously re-emitted, or transmitted unchanged, depending 
on the atomic state populations. A second beam split- 
ter 2 recomposes the two original field modes, and the 
photo-currents induced in detector Z)_|_ and Z3_ allow to 
measure the field phase shift. This measurement modi- 
fies the state vector of the two atomic clouds that become 



entangled, and at the same time it allows to probe the 
entanglement evolution. 




FIG. 2: Level structure of the atoms. The states |a) and \b) 
are stable states, \a) is coupled off-resonantly by the probe to 
the excited state |c). Here 7 is the spontaneous emission rate. 

Due to the spontaneous scattering, not all the photons 
will be revealed by the detectors and Therefore, 
in order to get a more realistic description of the pro- 
cess, it is important to evaluate the effects the scattered 
photons have on the atomic state vector evolution and on 
the entanglement dynamics. To this aim, we can imagine 
to detect the scattered photons by means of a spherical 
detector surrounding the two atomic clouds, and thus cal- 
culate how this detection, which corresponds to a spon- 
taneous scattering event, modifies the atomic wave func- 
tion (Fig. ^ . Notice that in order to preserve the purity 
of the atomic state we consider a spherical photodetec- 
tor revealing not only the presence of a scattered photon, 
but its (azimuthal) direction as well. 

The reset operator 

In this and in the next sub-section we will develop a 
formalism to account for the reset of the atomic sample 
wave function after the photons are detected. 

To this end, we first introduce a reset operator i?, 
whose action on the atomic state j^*) yields a new non- 
normalized state 1^'), affected by the spontaneous scat- 
tering photo-detection, that must be normalized to the 
final atomic state. The reset operator was first intro- 
duced in the context of the quantum jump approach to 
study the dynamical evolution of open quantum optical 
systems determined by continuous (gedankeiO measure- 
ments on the radiated field [H El 111 El iH IS- In 
deriving the analytic expression of the reset operator, we 
will consider the formalism introduced by Hegerfeldt E3 
and properly adapt it to our case study. 

Let us suppose that the photons are scattered one at a 
time by an atom belonging indifferently to one of the two 
samples. In the dipole and rotating wave approximations, 
the interaction Hamiltonian reads 

" fe,£ 
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where n — 1,2..., N goes over to the total number of 
atoms N, k and e are the photon wave vector and polar- 
ization component respectively, is the energy differ- 
ence between the atomic level |c) and \a), LUk is the field 
frequency, g^: ^ is the atom-photon interaction strength, 
and a/ is the annihilation operator for the field com- 
ponent interacting with the atoms. Henceforth we will 
consider h = c = I. Consequently ^ can be defined as 



(2) 



where dac is the transition dipole moment, the same for 
every atoms, e is the polarization vector and is the 
quantization volume. We are considering fields whose 
wavelength is larger than the spatial dimension of the 
two samples. Therefore, we can neglect not only the rel- 
ative positions of the atoms, but even that of the two 
atomic clouds, which are supposed to be placed in the 
origin of the coordinate system. As the photon interacts 
only with one atom at the time, we will consider this 
case in the development of the calculations. The gener- 
alization to the simultaneous interaction with N atoms 
is straightforward. 

During the interval At ^ the time evolution oper- 
ator, up to second order perturbation theory, is 



rt+At 

Ul{t + At,t) = dt'Hi{t') 
rt+At 



/t+At pt' 
dt' dt"Hi{t')Hi{t") 



(3) 



Let \i) = 



fc, 1^') be the initial state of the system 

photon-|-atom; after the interaction has taken place, it 

evolves in k' , e'^ l^fe/) = Ui{t + At, t) k, 1^-). 

Assuming that the photo-detection allows to discrim- 
inate only the direction of the photon wave vector, the 
correct final state is: 



|/) = ^|fc',e')(fc',e'|;7Ki + At,<) 



(4) 



From Eq. ^ it follows that I'i'k 

At,t) 
erator 



k',e' 



Ui{t 



k,e) 1^) and we can then introduce the reset op- 



k',e' 



Ui{t + At,t) 



k,e , , 



(5) 



acting on the atomic state and, once the photon is de- 
tected, yielding the updated atomic state vector l^'fc'). 



In our case, the term of first order in Eq. lO vanishes, 
and we have 



k',e' Ui{t + At,t) k,e)=6k,k'S. 



t+At ft' ,^ ^ V 

dt' dt"{k',e' Hi{t')Hi{t") k,ej. (6) 
Working out explicitly the second order term 



fc',e' 



Hj{t')Hi{t") 



k,e), we finally obtain 



k',e' 



Ui{t + At,t) 



t + At rt' 

dt' / (it"e'("'"-"''')*'e-*("''-"''')*" |c) (c| 
t Jt 



l-t+At nt' 



(7) 

As we are considering atoms initially in their stable 
states I a) and |6), the first integral in Eq. (0) acting on 
these states is always zero ((c|a) = (c|6) = 0). Thus, we 
are interesting only in the second integral. The calcula- 
tion is now similar to that of the transition amplitude for 
resonant scattering. Performing the substitution t — 
t + At = ^=>At = T, and adapting the method used in 
Ref. ^3 for the case of resonant scattering, we have 



k',e' Ui{t + At,t) k,e)^dk,k'Se.e' + 



2m- 



.9Le9y_. 



k',e' c-T 



5''{uk' - ujk) \a) {a\ 



(8) 



where A = cj/j — is the detuning between the in- 
cident field and the atomic transition energy, 7 is the 
spontaneous emission rate of the transition |c) — > \a), 
and the delta function 6'^, expressing conservation of en- 
ergy, reads S'^{uj — iv') — fj^ dr'e'^""" ^'^ = sin[(ijj — 
u;')]t/2tt{u;-uj'). 

li k' k and e' 7^ e, i.e. if the photon is scattered in a 
different direction than the incident one, Eq. © becomes 



k',e' Ui{t + At,t) k,e 



9k,e9p 

27ri— ^ r^(5'^(wfe/ - ujk) \a) {a\ . (9) 



Therefore, the non-normalized final state of the atom, 
after the photon is scattered and then detected, is 



-27ri 



9k,e9j:,,,, 



S^LUk' -tOk)\a){a\'<Sj) . (10) 
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From the above equation it is clear that, after the in- 
teraction with the photon, the atomic state is projected 
on the stable state \a), the same atomic state that was 
populated before the interaction. Of course, if initially 
the atom is in state \b), according to our assumptions it 
does not interact with the photon, and the latter is then 
transmitted unaffected. 

The probability that after the photo-detection the 
atomic state is |^fc'), is given by the square norm 
II I^Pfc') Ip. This coincides with the probability that a pho- 
ton is emitted in state fc',e'^ during a time interval r. 
Denoting ^,{^) = || 1*^-) |p, we have 



LUk'LOk 



A2 + ^ 



(2£oi3)2 

xSiLOk'-u;k)\\\a){a\^)\\'', (11) 

where we have written ^ according to Eq. ([SJl, and we 

have used the approximation liniT^oo [^'^{^k' ^ ^fe)]^ = 
TS{ujk' — LUk)/2'K. Because we are assuming that only the 
direction of the scattered photon is revealed, we have to 
trace over the polarization components and to sum on the 
energy spectrum. Thus, the probability that the photon 
is emitted along direction k' is 

t3 roc 



^ (27r)3 



cLv'cj'^I^,^^), (12) 



where \k\ = k = LOk- The integral in the above equation 
is easily evaluated, yielding 



totldac 



\dac 



;27r)24£2L3(A2 + :l) 



(am 



(13) 



In order to perform the sum over e', let us assume that 
the transition dipole moment is directed along the z axis, 
and that the wave vector of the incident photon points in 
the positive y's direction. Then we can write dac — dacZ, 



so that dr, 



since |ez| = 1, and dac ■ e' 



dac^'z- Moreover, making use of the standard relation [l^l 
Ee' Mac-e'P = dl^-dac-k'/ujl = ^^^(l-fcf), and taking 
k' = (sin 6* cos sin sin cos 61), where 6 and (p are the 
azimuth and polar angle respectively, we finally have 



cos 



(27r)24£§L3(A2 



2l\ 

4 / 



(am 



(14) 



It is now convenient for our purposes to express ly (5*) 
in terms of the interaction strength and of the sponta- 
neous emission rate. Actually, since 6^ — 1, from Eq. 
we can write gj: ^ — —dac^/^^k/'^SoL^ = g, and because 



— , ,3 



d^j,/37r£o, Eq. l(T^ becomes 
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Stt 
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(15) 



We can now write the explicit form of the reset oper- 
ator: 



(16) 



where = 37|g|2(l - cos^ 6')r/87r(A2 + 72/4). 

It is immediate to generalize this form to the case of 
an ensembles of atoms. We simply have 



N 



(17) 



and the probability that the photon is emitted along di- 
rection k' is 



4'W = II%I*> 



(18) 



Thus, after the photon has interacted with the atom, 
whose initial state is denoted by 1^*), and it has been 
detected, the updated atomic state, properly normalized, 
reads 



Ry \'^) 



(19) 



The effective time evolution operator 

Let us consider now the case in which the photon is 
emitted in the same state of the incident one, i.e., when 
k = k' and e = e'. Since for w « w' we have S'^{u! — u') « 
t/2-k, Eq. ^ can be rewritten as 



Ui(t + At,t) 



fc,e) =l-i- 



\g\' 



*2 



(20) 



If we take A ^ 7/2, to first order in 7/A2, we have 
1/(A - 17/2) w 1/A - i7/2A2. Adhering to this ap- 
proximation in Ea. l2()|) . we have 1 — i\g\'^T/ (A — 17/2) w 

1 - Z|5|V/A - |g|27r/2A2 W g-*|g|V/A-|5|^r/2A=^ 

Hence, the action of the time evolution operator on the 
state of the photon-atom system is simply to multiply 

this state by an exponential factor, i.e. 



,-i|ff|V/A-|g|V/2A2 



t). By extending this de- 
scription to the general case of A^ atoms, we can write 
the following effective non-Hermitian Hamiltonian 



Heff ^2^^a\ai\a)Aa\-i 



2A2 



a\di\a)^(a\, (21) 



where we have used 



fc, ey = a| |0). The physical mean- 
ing of the above expression is quite clear. The first 
term in the right-hand side is the phase-shift, propor- 
tional to the total population of state |a), caused by 
the off-resonant interaction between light and atoms. 
The second term represents the depletion of the field 
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mode a/ due to the spontaneous scattering. Thus, in 
the case the photon is ofF-resonantly transmitted, the 
atomic state evolution can be simply written as j^*') cx 

g(-i|s| = r/A-|g| V/2A2) E Ja)„(a| 

We can introduce a simplified description that eluci- 
dates the role of the probabilities for the alternative evo- 
lutions of the atomic state. Let us consider only one 
atom, whose initial state is \a) and, hence, will surely in- 
teract with the photon. If the photon is scattered in the 
same mode of the incident one, after an interval t, the 
atomic state is e"*'"' '^/^"Ifl 7t/2A |^^^ write 
the damping term as e"'^' jt/2A ~ ^ — \g\'^^T/2A'^ sa 
y/l — |gp7T/A2. Hence, we have (in the case of one 
atom) 



I*') 



(22) 



where || • || is the appropriate normalization factor. The 
factor in square in the above equation can be consistently 
interpreted as the probability amplitude that no sponta- 
neous scattering occurs. In fact, to first order in 
the probability that the photon is scattered in the k' 
direction by an atom in state \a) is: ^j^, = \\Iy{a)\\ — 
3|(7p7r(l — cos^ 9)/8ttA'^. By integrating over the whole 
solid angle, we recover the probability to have sponta- 
neous scattering: 



6 = 



\Iy{a)\\dn 



2-iT 



dip 



3 |5|V(1 



A2 



(23) 



Consequently, the above interpretation of the factor in 
square root in the right-hand side of Eq. H22|l follows. 

This result can be generalized by induction to the case 
of A'^ atoms. If Ua is the number of atoms populating 
state |a), using Eqs. (|17I18|I the probabihty of sponta- 
neous scattering can be written as = S.s'n^- Therefore 
the square root term becomes -^/l — Csn"^- Obviously, 
this approximation holds as long as ^s?^q ^ 1- 



III. CREATION OF ENTANGLEMENT BY 
PHOTO-DETECTION 

Let us introduce the atomic spin operators for the sta- 
ble states: = {\b)Jb\ - \a)Ja\)/2, j„+ = |6)„(a|/2, 
— \a)„{b\ /2. The reset operator lfT7|) and the eff'ec- 
tive Hamiltonian H21() become 



— ^ /^ 



(24) 



eff 



A 



alaj 



|gp7 

2A2 



Ed 



Jr. 



(25) 



Since we are concerned with ensembles of atoms in 
which each atom is initially prepared in the same state 
and interacts in an identical way with the surround- 
ing environment, the collective atomic state preserves 
the full permutation symmetry. It is therefore conve- 
nient to expand this collective state in the eigenstates 
of the effective collective angular momentum: \^) = 
'^M=-j ^m\J,M), where J — N/2 is the total angu- 
lar momentum, and M — {ni, — Ua) /2 is the eigenvalue 
of the operator Jz — J2n jnz- Collective raising and low- 
ering operators, and the corresponding Cartesian x- and 
y- components of the collective angular momentum, are 
defined as similar sums over all atoms in the sample. No- 
tice that if the atoms are prepared to populate only the 
two stable states \a) and \b), the total population of state 
\a) can be written as: J2n=i (1/2 - M) = {N/2 - M) = 
{ua + nb)/2 - [ub - na)/2 = Ua- 

The generalization of the atomic spin formalism to the 
case of two atomic samples with A^i and N2 atoms re- 
spectively, is straightforward. The initial disentangled 
state of the atomic samples is expanded in product- 
state wave functions of the two ensembles \'^) = 
Ea/i,m.^Mi,m. |Mi,M2), where |Mi,M2) = |Ji - 
Ni/2,Mi) (g) IJ2 = N2/2,M2), and we assume an ini- 
tial product state Ami,Mi — ■^MiAm2 before we let the 
system interact with the incident field. The total reset 
operator is simply the sum of the reset operators acting 
on a single ensemble of atoms as expressed by Eq. (|23J, 
Hence, we have 



A^i +iV2 



- (Ml + M2) 



\Mi,M2 



(26) 



where (iVi +iV2)/2- (Mi -t-M2) = na= AA(Mi+a/2) is the 
total population of the atomic state \a). 

After the photon has gone through the interferometric 
set-up and has been detected, the atomic state will be 
modified. In particular, if the state of the photon enter- 
ing the interferometer is = {\F)j^^ + |/)p^)/V2, 
where \F) , and denote the states of the pho- 
tons that follow respectively the upper or the lower path 
of the interferometer and |±)p^ = {\F)p^± \I)p^)/V2 
are the photon states revealed by the detectors _D+ and 
the initial state of the system photon-f atoms is 

l<^)ph+at = \'^)^i\F)ph + \I)ph)/V2. The analysis pre- 
sented in the previous sections provides us with the recipe 
to work out the evolution of the interacting component 
of \<l>)pii^af Therefore, by using Eq. (|26|l and the multi- 
atom generalization of Eq. (|22|l . after the photon has in- 
teracted, and before it is detected, we can formally write 
the photon+atoms state vector as 



/ ph+at 



E 

Ml, Ah 



A 



V2 



|i^)pJMi,Af2) 



6 



AMi,M; 



= E 

Ml, Ah 



1 + /i _ c p-«Ti-'rA/'(A/i+M2) 



l-./l-CsAA, 



where 15')^^ is the state of the photon scattered in the 

direction k' . Finally, to obtain the properly normalized 
state, we must sum over every possible spontaneously- 
scattered photon state. 

From the above equation is clear how the atomic state 
amplitudes change after the photo-detection. We must 
consider three cases: 

• If the photon is detected by I?+ or Z?_ , then the up- 
dating procedures for the atomic state amplitudes 
are respectively 



A 



■Mi,M2 



A 



I 1-^ 



(27) 



A 



A 



Ml, Ah 



All, Ah 



1 y^l ^sJ^^Mi+Ah) " 

and we get the non-normalized states 



(28) 



All, Ah 



1 ± ,/i - iMl 



(Afi + Ma) 



e ' A 



1 ,9 I 



(29) 



The probabilities that the photon is detected in Dj^ 
or £>_ are simply P± = || |^'±) |p, and the normal- 
ized atomic state is |*'^) = |*±) /|| |*±) ||. 



If the photon is spontaneously scattered, the atomic 
state amplitudes become 



A 



A 



Ml, Ah 



Ah -Ah 

71" 



^k'-^{Mi 



+Ah) 



(30) 



and the corresponding non-normalized state is 
1*^)= =^^V^AA(Ai,+A/2)|Mi,M2). (31) 



Ah, Ah 



In this case, the probability that the photon is scat- 



tered in the k' direction is Pf., 



and 



the normalized state reads l^'^) = l^'s) /|| |*s) |1- 
Let us notice that the probability that the pho- 
ton is scattered in any direction is simply Ps — 
J II |*s) W^d^}, and then P+ + P. + P, ^ 1. 

Form Eqs. H29|l and H31() . the entangled nature of the 
final atomic states is evident: because the amplitudes 
"evolved" into functions of the sum {Mi + M2), it is not 
possible to express the atomic state vector as a direct 
product of single atomic sample states. Every time a 
photon is detected, the updating procedure will be re- 
peated according the following simulation scheme: 

1. The probability Pg is determined and compared 
with a uniformly distributed random number 
(UDRN) n. 

2. If Ps < ri, then we calculate the conditional prob- 
ability that a photon is detected in Z)+, P"^ = 
P+/(l — Ps), and compare its value with an other 
UDRN rz: 

i) If Fiji > r2, then the atomic state amplitudes 

change according to Eq. lf77|) : 

ii) If < r2, then the atomic state amplitudes 

change according to Eq. 

3. If Ps > ri, we identify in which direction the pho- 
ton is scattered by determining the azimuthal angle 
0, whose density distribution function is given by 
2Tr^y/^s, and the atomic state amplitudes are up- 
dated using Eq. ^U^. 

4. Finally, we properly normalize the resulting atomic 
state vector and the process restarts. 

After Npii photons have been detected, if we denote 
with Ns the number of spontaneously scattered photons 
revealed by the surrounding detectors, and if N+ and 
are the photons detected by Z?+ and Z?_ respectively, 
such that Ns + -/V+ + = Nph, the atomic state vector 
will be 



I*') = 



1 



Ah, Ah 
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+ M2) 



X (Ml + M2)r_-{Mi + Ma) |Mi, M2) , 



where 



1 + h — f KT^ p-«^-r^(A/i+M2) 



and is the scattering factor corresponding to the n^- 
th detection along the different directions fc"s. 

IV. NUMERICAL SIMULATIONS 

In this section we will implement in a numerical simu- 
lation the detection model described in the preceding sec- 
tion. Following the measurement scheme introduced in 
Ref. ["llj , we will consider two atomic samples having the 
same number N of atoms and whose initial state is the 
eigenstate of the spin operator Jxi + Jx2, with eigenvalue 
J1 + J2 = 2 J = N. This means that all the atoms are ini- 
tially prepared in the superposition state (|a) -I- |6))/V2. 
The corresponding amplitudes in the basis eigenstates 
of the z-component of the collective angular momentum 
operators \J,Mi) are 



A 



Mi 




(2J)! 



(J + M,y.{j - M,y. 



(32) 



where i — 1,2. Expressed in terms of the number of 
atoms in state \a) in each sample, n^, the square of the 
amplitude ^ is I^a/.P ^ (l/2)^Ar!/(iV - <)!«)!. 
This means that in each ensemble the atoms are dis- 
tributed in state \a) and \b) according to a binomial dis- 
tribution with probability 1/2. 

In what follows, we will present the results of two 
kinds of numerical simulation where the entanglement of 
the two atomic samples is monitored. The first scheme 
corresponds to consecutive measurements of the angu- 
lar momentum operators Jzi + Jz2 and Jyi — Jy2, while 
the second corresponds to a simulated evolution in which 
continuous opposite rotations are applied to the atomic 
spin of the two samples as the photo-detection proceeds. 
Thereby we are considering a continuous exchange of the 
operators Jzi + Jz2 and Jyi — Jy2- We will compare the 
case in which the spontaneous scattering is neglected, i.e. 
5s = 0, with that in which the spontaneous emission rate 
7 is much smaller, but not negligible, than the detuning 
A. Specifically, we will take A = I5O7. 



Consecutive measurements of Jzi + Jz2 and Jyi — Jy2 

Starting from an initial state whose amplitude is given 
by the product of those of each samples expressed in 



Eq. H32() . we consider a series of photo-detections that 
project the state according to Eqs. and (pTT|l . 

As the photo-detection proceeds, the uncertainty in 
M12 = Ml M2, denoted with A2(Ji^ + J2z), goes to 
zero [T]| . i.e., after a large number of photons has been 
detected the state of the two ensembles approximates an 
eigenstate of Jzi + Jz2. This fact is true also if spon- 
taneous scattering is taken in account. Namely, we have 
seen that the detection of scattered photons measures the 
atomic state population given by (iVi -|- N2)/2 — (Mi -|- 
M2), just as the phase shift measurement carried out by 
detectors D+ and D_ does. If we take S = g^/A, we 
assume a phase angle St = 0.1, corresponding to a large 
(and hardly experimentally available) phase shift on the 
atomic state due to interaction with a freely propagating 
field. However, in our simulations a smaller and more re- 
alistic value of St just implies that more photons have to 
be detected to achieve the same reduction in the variance 
A2(Ji, + J2,). 

To quantify the pure-state entanglement between the 
two samples we will use the entropy of entanglement, de- 
fined by Bennett et al. to measure the entanglement 
of two systems in a pure state I^P): 

£ = -Tr{pi log2 pi) = -Tr{p2 \og^ 92), (33) 

where p\ = Tr2 |^) {^\ is the reduced density matrix of 
system 1, and analogously for p2. UN is the number 
of atoms in each sample, and we restrict ourselves to 
states which are symmetric under permutations inside 
the samples, the quantity £ takes values between zero 
for a product state, and \og2{N + 1) for a maximally 
entangled state of the two samples. 

Fig. O shows the evolution of entanglement between 
the two atomic samples obtained in four numerical simu- 
lations of the consecutive measurement of Jzi -I- Jz2 and 
Jyi — Jy2- For a number of detected photons Np^ < 2500, 
the state vector evolves according to the analysis exposed 
above, i.e., it goes approximately to an eigenstate of 
Jzi + Jz2- We then apply opposite rotations to the atomic 
samples and proceed with similar measurements as be- 
fore, which effectively measure Jyi — Jy2. After the rota- 
tions, the full permutation symmetry of the system is bro- 
ken, i.e., the total angular momentum = jj + + j2 
is not conserved. The final state of the two samples is 
still a simultaneous eigenstate of both Jf and J|, but 
with different values of the total angular momentum. 

Figs. EKa) and -(b) show the results of two simula- 
tions with spontaneous scattering. Fig. OKa) represents 
an example of the evolution of the entropy in which the 
final value {£Nph = 4.3165) is very close to that of the 
maximally entangled state (log2(20 + 1) = 4.3923). In 
Fig. |31-(b), we show a typical evolution, where the fi- 
nal value of the entanglement is well below the maxi- 
mum {£n^^ = 3.7494). Figs. EKc) and -(d) show the 
corresponding evolution of the entanglement when spon- 
taneous scattering is neglected. In Fig. |Sl-(c), the final 
value of the entropy is f jVp^ = 4.3194, while that reached 
in the simulation represented in Fig.|31-(d) is £Nph = 3.75. 
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(a) With spontaneous scattering (b) 




5000 



Number of detected photons 



FIG. 3: Entanglement of two samples each with = 20 
atoms. After detection of the first Nph ~ 2500 photons, the 
samples are rotated ±90 degrees in spin space, and the sub- 
sequent detection events lead, typically, to further entangle- 
ment. Figs, .(a) and (b) show results of two different simu- 
lation records for A = I5O7. Figs, (c) and (d) show results 
of two different simulation records for 7 = 0. In both cases, 
St = 0.1. 



The behavior of the two types of simulated evolu- 
tions are very similar, as confirmed in Fig. ^ where 
the average over 100 simulations for each type of in- 
teraction scheme is presented. Actually, in the first 
part of the simulation, for Nph < 2500, the two av- 
erage evolutions are practically the same, reaching the 
same average value of the entanglement for Nph = 2500 
{£n^^ = 2.7037). For Nph > 2500, after the rotations, 
despite the fact that the difference of the final average 
entropies (SNphl-y^o — SNphl-y^^o ~ 0-064) is smaller than 
the statistical fluctuations (AfjVphl7=o = 0.47), the aver- 
age value of the entanglement in simulations with spon- 
taneous scattering are constantly below those in which 
the spontaneous scattering is neglected. We will see that 
this feature is more evident in numerical simulations of 
measurements with continuous rotation of spin compo- 
nents. 

In the majority of the simulations, for both interaction 
schemes, the second round of measurements increases the 
entanglement. This fact can be qualitatively explained 
by noting that the approximate eigenstate of Jzi + Jz2, 
produced by the detection, has amplitudes on different 
Ml and M2 states with Mi -I- M2 fixed by the measured 
value, but the eigenspace is degenerate, and the ampli- 
tudes are simply proportional to the ones in the initial 
state. Due to the initial binomial distribution on Afi and 
M2, the distribution over, e.g.. Mi will therefore have a 
width of approximately ^fN . The reduced density matrix 
has the corresponding number of non-vanishing popula- 



tions, suggestive oi £ ^ \0g2vN = 0.51og2A^, half of 
the maximal value. This argument accounts for the first 
plateau reached in Fig. O 



FIG. 4: Average of the entanglement over 100 simulations of 
consecutive measurement of Jzi + Jz2 and Jyi — Jy2- Dashed 
line: average evolution when A = I5O7. Full line: average 
evolution when 7 = 0. Here 5t = 0.1 as in Fig.|3 



The measurement of the z-components causes a broad- 
ening of the distribution of the eigenstates of Jy, com- 
pared to the initial distribution, which was also binomial 
in that basis. The subsequent measurement of Jyi — Jyi 
will produce a state with M\y — Miy fixed by the mea- 
surement, but within the degenerate space of states with 
this fixed value the distribution on Mi^ of the reduced 



density matrix is broader than v N ^ and the entangle- 
ment is correspondingly larger. 



Measurements with continuous rotation of spin 
components 

The above arguments lead naturally to consider con- 
tinuous exchanges of the spin operator, in order to gener- 
ate states with higher entanglement. In this subsection, 
we present the results of simulations in which continuous 
opposite rotations are applied to the atomic spins of the 
two samples as the photo-detections proceed. If a is the 
variable rotation angle, what is effectively measured is 
the observable J a = cosQ!(J2i + Jzi) + sinQ!(Jj,i — Jyi)- 
If we denote with A^l^^^.j^ the wave function amplitude 
of the system after Nph photons have been detected, the 
updated state vector rotated by the angle ±a prior to 
the subsequent detection is 



I*; 



Nph X! X! '^MiM[ ^'^ 



Mi,M2 



Ji, M- ) is the ro- 



where D'^^^^^, {a) = (M„ J,| e-^''-(") 

tation matrix element for sample i {i — 1,2). The wave 
function updating algorithm has the same structure as 



9 



the one described in section Ijllljl . In a realistic experi- 
mental situation, one would detect photons according to 
a Poisson process. With a constant rotation frequency in- 
duced, e.g., by applying opposite DC magnetic fields onto 
the atoms, this would lead to small rotation angles with 
an exponential distribution law. For simplicity, however, 
we apply the same small rotation angle a between each 
detection event. 

With spontaneous scattering 
4 




a. 
o 




5000 

Wittiout spontaneous scattering 

4 



E 



2500 5000 

2500 5000 "0 2500 

Number of detected photons 




FIG. 5: Entanglement of atomic samples with 20 atoms each. 
The spins are rotated in opposite directions around the a:-axis 
by the angle a — n/5 after each detection event. Insets (a- 
d) show the evolution of the overlap | {5'o| ^')iVj,^P between 
the state of the samples and the maximally entangled state. 
Figs, (a) and (b) present results of two different simulations in 
the case A — I5O7. Figs, (c) and (d) show analogue results for 
simulations in which 7 = (no spontaneous scattering). In (a) 
and (c) ^ converges toward I'I'o) and further detections 

have no effect on the state vector. In (b) and (d) the state 
gradually loses its component along I'I'o), and it subsequently 
evolves in the orthogonal subspace of |^o). Here, as in the 
previous simulations, St — 0.1. 



In Fig.|Slwe show four numerical simulations of the en- 
tanglement evolution, two for each interaction schemes. 
The rotation angle between subsequent photo-detection 
events is a = tt/S. In the upper part of the panel, we 
provide results of two numerical simulations with sponta- 
neous scattering. Fig.|31-(a) shows the evolution in which 
the maximal value of the entanglement is reached, as is 
confirmed by the inset of the figure showing the evolution 
of the overlap between the state vector of the two atomic 
samples and the maximally entangled state 



(34) 



In Fig. ]^{h) we have a typical case in which the final 
state is orthogonal to |^o)- The lower figures has been 
obtained from simulations in which the spontaneous scat- 
tering has been neglected. In Fig. \^{c) the evolution 



leading to the maximally entangled state is shown. A 
typical result is shown in Fig.|Sl-(d). 

Note that 1^*0) is the only joint eigenstate of the opera- 
tors Jjx — J2xi Jiy — J2y and Jiz + J2z with null eigenvalue 
hllm. When we simulate the detection of phase shifts 
proportional to Jiz+J2z and Jiy — J2y, or combinations of 
these operators, there is a nonvanishing probability that 
the state vector is gradually projected onto |^'o), and this 
state is unaffected by all further photo-detection events, 
see Figs.EKa) and -(c). 

If the state of the atomic samples has not collapsed into 
1 5*0) after a large number of photons have been detected, 
it instead gradually becomes orthogonal to that state, 
and the maximally entangled state is never attained, see 
the insets in Figs.EKb) and -(d). 




1000 2000 3000 4000 

Number of detected photons 



5000 



FIG. 6; Average of the Entanglement over 100 simulations of 
measurements with continuous rotations of the spin compo- 
nents. Dashed line: average evolution when A = I5O7. Full 
line: average evolution when 7 = 0. Here 5t = 0.1, as in the 
previous figures. 

In these simulations the qualitative behavior of the 
two types of interaction schemes is again similar. As ex- 
pected, there is a faster increase of the entanglement and 
a more efficient evolution toward a maximally entangled 
state, compared with the simulations of consecutive mea- 
surements of Jiz + J2Z and Jiy — J2y When spontaneous 
scattering is considered, the evolution toward a stable 
regime appears to be slower, as is evident in Fig. |^ where 
an average over 100 simulations is shown. Moreover, in 
this case it appears clearly that the average values of the 
entanglement for evolutions with spontaneous scattering 
is constantly smaller than that obtained from evolutions 
without spontaneous scattering. Even though this effect 
is very weak, it seems to become evident when the rota- 
tional symmetry of the two systems is broken, as seen in 
the previous section, and it is suggestive of less correlated 
atomic samples due to photo-detections of spontaneous 
scattered photons. Namely, when spontaneous scattering 
is considered, we have a wide spectrum of possible states 
in which the system can be found. Hence, in detecting 
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all the scattered photons, we gain information from a 
more disordered system in comparison with the case in 
which spontaneous scattering is neglected. In the latter 
case there are only two possible states Tl^ and less igno- 
rance (the entropy is smaller) about the system. Yet, the 
good news is that the effect of spontaneous scattering is 
not dramatically destructive: our analysis provides indi- 
cation that a QND interferometric scheme to effectively 
entangle atomic samples should work very well also when 
this unavoidable physical process is taken into account. 

It is worth noting that even when the final state does 
not coincide with |4'o), we observe an entanglement that 
is almost constant as the photo-detection goes on. This 
suggests the presence of families of states, orthogonal to 
|4'o), with relatively well-defined values of the operators 
Jiz + J2z and Jiy — J2y, which is allowed by Heisen- 
berg's uncertainty relation as long as the expectation 
value of the operator Ji^ — J2x is small This inter- 
esting threshold behavior is unexpected, and at present 
not completely understood. 

V. CONCLUSIONS AND OUTLOOK 

In this paper we have analyzed the effects of the spon- 
taneous scattering on the entanglement evolution of two 
atomic samples correlated by QND measurements of the 
atomic state population by means of optical phase shifts. 

We generalized the interferometric set-up introduced 
in Ref. to measure the field phase shift, by adding 
a photo-detector surrounding the two atomic clouds to 
allow detection of the scattered photons. In this way 



we were able to study the effects of spontaneous scat- 
tering, neglected in the ideal treatment We have 
introduced a formalism of conditional quantum dynam- 
ics in terms of reset operators that, acting on the initial 
atomic wave functions, yield the new state vector modi- 
fied according to the result of the photo-dctcction of the 
spontaneously scattered photons. We have then deter- 
mined the associated effective Hamiltonian that provides 
the time evolution of the system in the instances that the 
photons reach the detector measuring the phase shifts. 
We have finally implemented the formalism in numeri- 
cal simulations, by means of which the evolution of the 
entanglement of the two atomic samples has been moni- 
tored. We have compared the results of the simulations 
in the presence of spontaneous scattering with the ideal 
evolution that neglects this effect. Under the assump- 
tions and approximations adopted, we have shown that 
the effects of the spontaneous scattering are rather small, 
and thus that the QND interferometric scheme should be 
robust in realistic situations. 

The model dynamics introduced in the present work 
should be useful also for further studies on the dynam- 
ics of entanglement of atomic samples. In particular, it 
would be very interesting to combine the entangling pro- 
tocol presented here with some feedback scheme [23, , 
in order to drive the system with very high probability 
toward the maximally entangled state. Moreover, the 
conditional dynamics could be used to investigate the 
entanglement properties of the state space in which the 
atomic system state vector is projected by the measure- 
ment, with the scope of exploiting these properties to 
realize more efficient entangling protocols. 
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